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S joberg , Wi 1 1 iam  F..M.S.,  Purdue  University,  May  1986. 
Modification  of  Parameterized  Latent  Heat  Release 
Estimates  Using  Unenhanced  and  Enhanced  Satellite  Imagery. 
Major  Professor:  Dr.  Phillip  J.  Smith 

This  study  investigates  the  use  of  unenhanced  and 
enhanced  satellite  imagery  to  modify  parameterized  latent 
heat  release  for  a  developing  cyclone  over  the  central 
United  States  during  10-11  April  1979.  Parameterized 
stable  and  convective  latent  heat  release  estimates  are 
evaluated  initially  by  a  qualitative  comparison  with 
National  Weather  Service  radar  summaries:  gridded,  time- 
averaged  observed  precipitation  fields;  and  satellite- 
derived  cloud  cover.  This  evaluation  is  then  extended  to 
the  latent  heating  values  obtained  after  applying  the 
three  modification  schemes.  Finally,  the  latent  heating 
values  are  evaluated  quantitatively  using  the  gridded, 
observed  precipitation  values.  This  evaluation  includes 
comparisons  of  area-averaged  totals  and  correlation 
coefficients  for  the  total  SESAME  grid.  In  addition,  the 
total  number  of  grid  points  at  which  observed 
precipitation  and  latent  heating  values  occurred  over  both 
the  total  SESAME  grid  and  three  subregions  representing 


different  air  mass  regions  are  compared. 

Results  indicate  that  the  satellite  schemes  yield 
improved  latent  heat  release  estimates.  This  is  seen  in 
the  improved  areal  coverage  of  modified  latent  heat 
release  fields  when  compared  with  radar  and  observed 
precipitation  summaries.  In  addition,  a  21Z-45Z 
improvement  in  the  average  correlation  between  calculated 
LHR  and  observed  precipitation  values  is  obtained. 
Furthermore,  the  number  of  grid  points  at  which  LHR  was 
calculated  matches  the  observed  precipitation  grid  points 
more  closely  after  the  modifications  were  applied.  When 
the  three  modifications  are  compared,  it  is  seen  that 
there  is  considerable  support  for  identifying  the  second 
modification  as  yielding  the  best  match  between  latent 
heat  release  values  and  the  comparative  data  forms.  This 
modification  uses  a  235°K  threshold  to  differentiate 
between  convective  and  stable  regions  of  a  cloud  mass.  It 
was  also  found  that  the  impact  of  these  modifications 
varies  depending  on  the  accuracy  of  the  initial  latent 
heat  release  calculations. 
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1.  INTRODUCTION 


Relationships  between  atmospheric  circulation  systems 
and  latent  heat  release  occurring  during  cloud/ 
precipitation  formation  have  been  a  matter  of  major 
concern  in  meteorology  for  a  number  of  years  (e.g. , 
Aubert ,  1957;  Danard,  1964.1966;  Tracton , 1973 ;  Vincent  et 
a  1 .  . 1 977 ;  Robertson  and  Smith, 1983;  Gyakum,  1983a, b;  Smith 
et  a  1 .  . 1984) .  A  typical  diagnosis  of  latent  heat  release 
involves  partitioning  the  quantity  into  two  components, 
stable(grid-scale)  and  convective  (subgrid-scale).  The 
stable  component  is  generally  determined  as  a  function  of 
the  grid-scale  vertical  motion  and  the  vertical  gradient 
of  water  vapor.  Methods  of  quantifying  the  convective 
component  have  concentrated  on  two  approaches.  The  first 
approach  focuses  on  the  microphysical  processes  of 
individual  cumulus  clouds  or  cloud  ensembles.  Schemes 
utilizing  this  approach  have  been  advanced  by 
Kess 1 er ( 1969) ,  Simpson  and  Wiggert ( 1971 ) ,  Takeda( 1971 ) , 
and  Orville  and  Kopp(1977).  The  second  approach  centers 
on  the  parameterization  of  subgr id-sca 1 e ( cumu 1  us ) 
processes  using  available  grid-scale  information  .e.g., 
parameterizations  by  Kuo( 1965 , 1974) ,  Ooyama( 1971 ) ,  Arakawa 
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and  Chao(1973),  and  Krishnamurti  and  Moxim(1971).  A 
general  discussion  and  comparison  of  the  various  cumulus 
parameterization  schemes  can  be  found  in  Cese 1  ski ( 1973) 
and  Krishnamurti  et  a  1 .  (1980).  The  scheme  used  in  this 
study  was  introduced  by  Kuo ( 1965 , 1974) ,  and  later  modified 
by  Edmon  and  Vincent ( 1976)  and  Lin  and  Smith(1979). 

Although  much  has  been  done  to  produce  realistic  latent 
heat  estimates,  diagnostic  studies  to  date  reveal 
deficiencies  in  both  the  magnitudes  and  the  distribution 
of  latent  heating  values.  Since  latent  heating  is  closely 
tied  to  the  occurrence  of  clouds  and  precipitation, 
information  about  the  latter  provides  independent  data 
which  can  be  used  to  improve  previously  obtained  estimates 
of  latent  heating.  In  this  context,  earth  orbiting 
satellites  are  particularly  useful  because  they  provide 
platforms  for  estimating  cloud  conditions  over  large 
areas.  Further,  as  suggested  by  Barrett( 1970 , 1971 ) , 
satellite  data  can  be  used  to  estimate  rainfall  amounts. 
Thus,  it  should  be  possible  to  use  satellite  information 
to  improve  latent  heat  release  estimates. 

There  has  been  considerable  work  done  on  the  use  of 
satellite  data  to  estimate  rainfall  amounts.  Barrett  and 
Martin(1981)  divided  these  methods  into  three  categories: 
1)  cloud-indexing  methods,  (Barrett,  1970,1971);  2)  life- 
history  methods,  (Griffith  et  a  1 ■  ,  1976,  1978;  Woodley  et 
a  1 ■  .  1980;  Scofield  and  Oliver,  1977);  and  3)  passive- 
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microwave  methods,  (Wilheit,  1972;  Smith  and  Kidder, 
1978).  Each  of  these  methods  acknowledges  the  importance 
of  accurate  satellite  cloud  estimations. 

In  the  present  study  the  utilization  of  satellite  data 
is  taken  to  its  next  logical  step,  that  of  combining 
direct  estimates  of  latent  heating  with  satellite  data  in 
order  to  compensate  for  the  inherent  deficiencies 
contained  in  each.  Latent  heat  release  calculations  are 
easily  made  from  the  synoptic-scale  data  network,  yet  can 
suffer  from  inadequate  data  coverage,  erroneous  data,  and/ 
or  computational  model  deficiencies.  Satellite  data  give 
spatially  complete  and  detailed  views  of  cloud 

distributions,  but  do  little  to  account  for  the  large 
changes  in  rainfall  found  in  small  regions  of  convective 
I  clouds,  precipitation  variations  that  can  occur  within 

similarly  appearing  cloud  fields,  or  the  vertical 

distribution  of  the  resulting  latent  heat  release. 

The  work  presented  here  is  an  extension  of  the  simple 
technique  described  by  Lin  and  Smith(1985),  in  which 
calculated  latent  heating  is  modified  using  satellite- 
derived  cloud  information.  Using  enhanced  infrared 
satellite  imagery,  the  Lin  and  Smith  procedure  is  refined 
in  two  ways,  one  to  permit  the  discrimination  between 
stable  and  convective  cloud  forms  and  the  other  to  allow  a 
determination  of  the  heights  of  satellite-derived  cloud 
tops  . 
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2.  DATA  AND  COMPUTATIONS 


2. 1  Rawinsonde  data 

The  period  of  investigation  consists  of  seven  synoptic 
times,  3h  apart,  from  1200GMT  10  April  through  0600GMT  11 
April  1979.  These  times  were  chosen  to  correspond  to  a 
much  publicized  outbreak  of  severe  weather  in  the  central 
United  States.  An  extratrop i ca 1  cyclone  with  significant 
convective  activity  and  associated  latent  heat  release  was 
observed  by  the  special  data  network  initiated  for  the 
AVE-SESAME  (Atmospheric  Variability  Experiment-Severe 
Environmental  Storms  and  Mesoscale  Experiment)  field 
program.  The  network  consisted  of  23  National  Weather 
Service  (NWS)  and  16  supplementary  rawinsonde  sites 
(Fig.2.1),  with  observations  taken  at  3h  intervals  over  a 
24h  period  from  1200GMT  10  April  -  1200GMT  11  April  1979, 
a  period  commonly  known  as  SESAME  I.  Further  details  of 
the  AVE-SESAME  '79  experiment  are  given  by  Alberty  et  a  1 , 
(1979),  Hill  et  al  .  (1979),  and  Gerhard  et  a  1 .  (1979). 

The  observational  data  set  included  geopotential 
height,  temperature,  relative  humidity,  horizontal  wind 
components,  and  surface  pressure.  These  data  were  then 
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objectively  analyzed  using  the  Barnes  (1964,1973)  scheme 
on  a  12  X  18  1  at i tude/ 1 ong i tude  grid  with  one  degree 
latitude/ 1 ongitude  resolution  (Fig. 2. 2).  All  data  were 
interpolated  to  50  mb  levels  extending  from  1000  mb  to  100 
mb.  Including  the  surface,  this  gave  a  total  of  21  levels 
of  data.  The  gridded  data  set  was  provided  to  the  author 
by  Dr.  F.  R.  Robertson  of  the  NASA  Marshall  Space  Flight 
Center . 

Vertical  motions  were  estimated  from  horizontal  wind 
components  by  the  kinematic  method,  which  utilizes  the 
continuity  equation  in  pressure  coordinates  (O'Brien, 
1970).  Because  there  are  systematic  errors  inherent  in 
the  evaluation  of  divergence  fields,  unrealistically  large 
estimates  of  w  are  often  produced  at  the  upper  levels. 
Such  bias  errors  were  suppressed  by  the  quadratic 
adjustment  scheme  suggested  by  O’Brien  (1970)  and  applied 
by  Smith  ( 1971 ) . 


2 . 2  Precipitation  data 

Observed  precipitation  rates,  used  to  evaluate  the 
vertically  integrated  values  of  parameterized  latent  heat 
release,  were  obtained  from  the  Hourly  Precipitation  Data 
booklets  published  by  the  National  Climatic  Data  Center, 


NOAA .  Precipitation  data  were  smoothed  to  the  grid  by 
averaging  reported  pr ec i p i ta t i on  quantities  over  each  one 
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degree  grid  square  for  3h  periods,  i.e.,  lh  before  to  2h 
after  the  seven  synoptic  times  under  consideration.  These 
data  were  provided  to  the  author  by  Dr.  D.  G.  Vincent, 
Department  of  Geosciences,  Purdue  University. 

2 . 3  Calculation  of  latent  heat  release 

Convective  latent  heat  release  (CLHR)  was  computed 
using  Kuo's  (1965,1974)  parameterization  scheme  as 
modified  by  Edmon  and  V-ncent  (1976)  and  Lin  and  Smith 
(1979).  Stable  latent  heat  release  (SLHR)  was  evaluated 
from  a  scheme  used  by  Krishnamurti  and  Moxim  (1971).  Both 
heating  components  were  allowed  to  occur  simultaneously  in 
the  same  region. 

In  the  Kuo  scheme  the  CLHR  is  parameterized  by: 

g( 1— b)L  Mt[  Tc(p)  -  T(p))  0/T 

CLHR*  -  (1) 

(  pb  -  pt)  <  Tc  -  T  > 

for  those  levels  with  Tc  >  T,  where 
1  P«  iq 

bMt  -  "  j*  —  dp  . 

g  pu  bt 

1  Pb 

<  Tc  -  T  >  *  -  f  (  Tc  -  T  )  dp  . 

Pb  "  Pt  Pt 


.  V 


Mt  *  the  rate  of  moisture  convergence  from  the  surface  to 
the  cloud  top,  including  horizontal  convergence  of 
specific  humidity  and  evaporation  of  water  from  the 
lower  surface; 

T  *  the  environmental  temperature; 

Tc  =  the  temperature  within  the  cloud,  assumed  to  be  the 
temperature  of  the  moist  adiabat  extending  from  the 
cloud  base; 

0  “  the  potential  temperature; 

L  =  the  latent  heat  of  condensation, 

*  2. 501x10s  JKg" 1  for  T  £  -20°  C, 

=  2.835x10s  JKg"1  for  T  <  -20°  C; 
pb  =  the  pressure  at  the  cloud  base,  the  base  of  the  first 
conditionally  unstable  layer; 

pt  =  pressure  at  the  cloud  top,  the  level  at  which  the 
moist  adiabat  from  the  cloud  base  crosses  the 
environmental  sounding;  and 

b  =  the  fraction  of  Mt  stored  in  the  air  to  increase  the 
specific  humidity  of  the  environment. 

The  Kuo  scheme  contains  the  following  modifications 
(Edmon  and  Vincent , 1976 ;  Lin  and  Smith, 1979).  A  latent 
heating  value  was  calculated  only  if  1)  the  cloud  was  at 
least  100  mb  in  depth,  2)  the  mean  vertical  motion  within 
the  cloud  layer  was  upward,  and  3)  the  mean  relative 
humidity  between  the  surface  and  500  mb  (RH.V)  was  greater 
than  0.4.  The  partitioning  parameter  b  was  evaluated  as  a 
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function  of  RH.V  in  a  fashion  similar  to  Anthes(1977)  and 
Donner  and  Kuo(1982).  If  the  value  of  RH.V  was  less  than 
or  equal  to  0.4,  then  b  was  assumed  to  be  one  (i.e.  .  no 
convective  latent  heat  release).  If  RH,V  was  greater  than 
or  equal  to  0.8,  then  b  was  assumed  to  be  zero.  This 
upper  bound  corresponds  with  the  relative  humidity 
constraint  imposed  on  the  stable  precipitation.  If  RH.V 
was  between  0.4  and  0.8,  b  was  assumed  to  be  l-RHav. 

SLHR  was  calculated  from  the  equation, 

SLHR-  -Lw  -  ,  (2) 

dp 

where  u  is  the  vertical  velocity  in  pressure  coordinates 
and  qs  is  the  saturated  specific  humidity.  Stable  latent 
heat  release  was  specified  under  the  following 
limitations;  1)  grid-scale  upward  vertical  motion,  2) 
vertical  moisture  convergence  for  saturated  air  at 
environmental  temperatures,  3)  relative  humidity  greater 
than  or  equal  to  80%  (Krishnamurti  and  Moxim,  1971;  Edmon 
and  Vincent,  1976). 

Henceforth,  the  sum  of  CLHR  and  SLHR  will  be  identified 
as  total  latent  heat  release  (TLHR) .  Results  obtained 
from  these  calculations  are  referred  to  as  the  "basic" 
latent  heat  release  (BLHR). 


2.4  Modification  of  latent  heat  calculations 
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Satellite  images  were  obtained  for  the  seven  time 
periods  under  consideration.  The  modification  schemes 
were  based  on  both  unenhanced  and  enhanced  infrared 
satellite  images.  Only  IR  data  were  used  in  order  to 
standardize  the  modification  schemes  for  use  during  both 
day  and  night.  Using  the  basic  analysis  grid  described 
earlier,  each  grid  square  was  assigned  a  sate  1 1 ite— derived 
cloud  fraction  based  on  the  percent  cloud  cover  estimated 
for  that  square.  The  satellite  cloud  fraction  value  was 
used  to  modify  the  basic  latent  heat  estimates.  Earlier 
work  done  by  Lin  and  Smith(1985)  describes  the  use  of 
unenhanced  satellite  data  to  modify  latent  heat  release 
calculations.  This  study  applies  their  technique  to  the 
SESAME  I  case  and  extends  their  work  by  including  enhanced 
sate  1 1 ite  data . 


2.4.1  Description  of  satellite  cloud  data 

All  satellite  images  were  derived  from  the  SMS-1 
satellite  observing  system  and  were  provided  to  the  author 
by  the  NASA  Marshall  Space  Flight  Center.  During  the 
SESAME  I  period  the  SMS-1  was  in  a  geosynchronous  orbit  at 
an  approximate  altitude  of  35800km  over  the  equator  and 
was  located  at  a  fixed  longitude  of  75°W.  The  SMS-1  was 
equipped  with  a  visible  and  an  infrared  spin  scan 
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radiometer (VISSR)  which  scans  the  full  disc  of  the  earth 
in  18.2  minutes,  viewing  in  the  visible  spectrum  (0.55  to 
0.754m)  and  in  the  infrared  spectrum(  10 . 5  to  12.64m). 
Visible  and  infrared  data  resolutions  are  approximately 
4km  at  the  satellite  subpoint.  Subsequent  sections  will 
discuss  the  logic  behind  the  specific  enhancements  used  in 
this  study  and  how  these  enhancements  were  incorporated 
into  the  latent  heat  release  modification  schemes.  More 
details  regarding  the  SMS-1  satellite  and  its 
instrumentation  can  be  found  in  National  Space  Science 
Data  Center (NSSDC)  publication:  "Data  Catalog  Series  for 
Space  Science  and  Applications  Flight  Missions:  Volume 
4a,”  Ng  and  Shue  editors. 

2.4.2  Lin  and  Smith  modification 

Using  unenhanced  satellite  imagery,  for  each  grid 
square  a  satellite-derived  cloud  fraction  value  of  0,  0.5, 
or  1  was  assigned.  A  value  of  1  was  given  to  those 
squares  with  over  80%  of  its  area  covered  with  cloud. 
Squares  with  less  than  20%  of  its  area  covered  with  cloud 
were  assigned  a  value  of  0.  For  squares  with  satellite 
cloud  coverage  between  20%  and  80%,  a  value  of  0.5  was 
assigned.  These  satellite  cloud  fractions(a)  were  used  to 
modify  the  original  total  latent  heat  release  estimates  at 
every  grid  point  and  level  using  the  following  procedure, 
tf  at  a  gridpoint(o)  the  original  latent  heat 


oa 1 cu 1  at  ion (Q0 )  was  non-zero,  a  modified  latent  heat 
est imate (Q0 * )  was  determined  by 

Q0*  -  «Q0.  (3) 

If  Q0  was  determined  to  be  zero,  and  a«0 ,  then 

Q0*  -  0.  (4) 

If  Q0  was  determined  to  be  zero  but  a*0 . 5  or  1.0,  then  the 
modified  value  was  determined  by  inverse  distance  weighted 
interpo 1 ati on 

1 

Z -  Q, 

d£ 

Qo*  "  -  .  (5) 

1 

I - 

di 

where  ds  is  the  distance  from  grid  point  o  to  a 
surrounding  grid  point  i.  For  grid  points  located  at 
least  two  grid  squares  from  the  computational  boundary, 
the  interpolation  was  done  using  points  two  grid  squares 
in  all  directions,  a  total  of  24  interpolating  points. 
Points  located  one  grid  square  from  the  boundary  used  the 
immediately  surrounding  grid  points,  a  total  of  8 
interpolating  points.  The  points  located  on  the  boundary 
used  the  3  adjacent  computational  grid  points. 

The  modification  scheme  was  applied  to  the  sum  of  the 
convective  and  stable  latent  heat  release  (TLHR)  and  will 
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be  referred  to  in  later  chapters  as  modified  latent  heat 
release  one  (MLHR1). 


2.4.3  235°K  enhancement  modification 

The  next  step  in  this  inquiry  was  an  attempt  to  isolate 
the  convective  region  of  satellite  cloud  masses  by 
enhancing  the  satellite  images  using  235°K  as  a  threshold 
temperature.  A  comparison  done  by  Arkin(1979)  of  6h 
rainfall  accumulations  over  the  GARP  Atlantic  Tropical 
Experiment  (GATE)  B-scale  array  and  fractional  cloud 
coverage  below  various  IR  threshold  temperatures,  showed 
that  the  maximum  correlation  of  these  two  factors  occurred 
at  235°K.  It  was  found  that  despite  GATE's  phase-to-phase 
variablity  in  factors  such  as  average  rainfall  rate 
(Hudlow  et  a  1 .  ,  1979)  and  the  mean  surface  wind  field 
(Seguin  et  a  1 .  .  1978),  the  relationship  between  rainfall 
and  the  fractional  IR  cloud  coverage  defined  by  this 
particular  threshold  temperature  remained  nearly  the  same 
throughout  GATE. 

Each  square  in  the  grid  was  divided  into  regions  with 
cloud  tops  colder  than  235°K,  warmer  than  235°K,  or  cloud 
free.  Modifications  to  the  basic  latent  heating  values 
were  made  by  assuming  that  tops  235°K  or  colder  contained 
only  convective  latent  heat  release  and  tops  warmer  than 
235°K  contained  only  stable  latent  heat  release.  The 


fraction  of  each  cloud  type  was  estimated  for  each  grid 
square  and  was  utilized  to  modify  the  basic  latent  heating 
values  in  the  following  way.  In  grid  squares  where  there 
was  less  than  20%  of  both  types  of  cloud,  it  was  assumed 
that  no  latent  heat  release  had  occurred  and  BLHR  values 
were  set  to  zero.  If  either  type  of  cloud  filled  greater 
than  80%  of  the  square  it  was  assumed  that  only  SLHR 
(CLHR)  occurred  if  the  cloud  top  was  warmer(equal  to  or 
colder)  than  235°K.  In  this  case  the  non— occurring  basic 
latent  heating  component  was  set  to  zero.  For  grid 
squares  having  some  tops  warmer  than  this  threshold  while 
other  tops  were  colder,  the  latent  heat  release  values 
were  adjusted  in  the  following  manner.  If  the  square  was 
covered  by  20%  to  80%  of  one  type,  while  the  percentage  of 
cloud  of  the  other  type  was  less  than  20%,  the  modified 
TLHR  value  was  assumed  to  be  half  of  the  basic  SLHR(CLHR) 
value  if  the  predominant  cloud  top  was  warmer(equal  to  or 
colder)  than  235°K.  If  there  was  between  20%  and  80%  of 
both  types  of  cloud,  the  modified  TLHR  was  assumed  to  be 
half  of  the  BLHR  value. 

Finally.  if  there  was  greater  than  20%  of  any  type  of 
cloud  and  there  was  no  latent  heat  release  calculated  for 
that  grid  square,  the  initial  zero  values  of  both 
components  were  altered  using  the  inverse  distance 
weighted  interpolation  procedure  summarized  in  the 
previous  section  and  then  the  SLHR  and  CLHR  values  were 


modified  as  described  above.  The  resulting  latent  heat 
release  values  from  this  modification  will  be  referred  to 
in  the  results  chapters  as  modified  latent  heat  release 
two (MLHR2) . 

2.4.4  3-layer  enhancement  modification 

The  previous  two  schemes  made  no  attempt  to 
discriminate  the  various  heights  of  observed  clouds;  both 
assumed  that  any  observed  cloud  extended  throughout  the 
1000  to  100  mb  layer.  Thus,  a  third  modification  was 
devised  which  differentiated  between  low,  middle,  and  high 
clouds.  This  differentiation  was  accomplished  by 
selecting  three  different  temperature  threshholds  for 
enhancement  of  the  infrared  photographs.  By  focusing  on 
the  vertical  extent  of  the  cloud  it  was  hoped  that  the 
latent  heat  release  calculations'  could  be  modified  to 
ensure  that  latent  heat  release  was  confined  to  levels 
where  cloud  was  present. 

The  definitions  of  low,  middle,  and  high  clouds  were 
adopted  from  the  Weather  Almanac  (1977),  which  specifies 
the  height  ranges  shown  in  Table  2.1.  Combining  these 
definitions  with  the  equivalent  blackbody  temperatures ( °K ) 
for  typical  cloud  levels  identified  by  Barrett ( 1972) , 
summarized  in  Table  2.2,  a  range  of  temperature  values  was 
identified  which  made  it  possible  to  separate  the  various 


clouds  by  level.  To  test  the  applicability  of  these 
temperature  ranges  to  the  current  data  set.  six  grid 
points  were  chosen  and  vertical  temperature/dewpoint 
soundings  were  anal  2ed  for  two  times,  1200GMT  10  April 
1979  and  0300GMT  11  April  1979.  These  points,  identified 
by  an  X  on  Fig  2.2,  were  chosen  in  order  that  each  of  the 
air  masses  involved  in  this  synoptic  disturbance  would  be 
represented.  Table  2.3  shows  the  temperature  ranges 
observed  at  these  points  for  each  of  the  height  categories 
specified  in  Table  2.1.  These  results  show  a  strong 
resemblance  to  the  temperature  values  observed  by 


Barrett( 1972) 


Thus,  the  Barrett  values  represent 


reasonable  threshold  temperatures  for  the  enhancement  of 
the  IR  satellite  images  for  this  data  set  The  temperature 
ranges  and  corresponding  color  code  for  the  enhanced  IR 
images  for  each  cloud  category  are  summarized  in  Table 
2.4.  The  significant  warming  at  the  surface  over  the  grid 
late  in  the  study  period  was  captured  by  a  reduction  in 
the  maximum  temperature  threshold  used  for  designating  low 
c 1 oud . 


To  provide  data  for  the  latent  heat  modification  a  grid 
was  overlayed  on  the  color-coded  satellite  images  and  the 
fraction  of  each  cloud  category  was  estimated  for  each 
grid  square.  The  predominant  cloud  level  determined  the 
final  modification  of  latent  heat  release  values.  For 


example,  if  a  grid  square  was  covered  with  50%  high  cloud 


Table  2.1.  Cloud  level  definitions  for  MLHR3. 


Table 


Table 


Tab  1  e 


CLOUD  LEVEL 
LOW 
MIDDLE 
HIGH 


HEIGHT  OF  BASES! KM) 
SFC  to  2 
2  to  7 
7  to  14 


2.2.  Equivalent  blaekbody  temperature ( °K)  for  MLHR3 
cloud  level  determinations  (from  Barrett 
, 1972) . 


CLOUD  LEVEL 


LU 

MIDDLE 

HIGH 


EQUIVALENT  BLACKBODY  TEMPftO 
200K  -  300K 
25DK  -  230K 
250K 


2.3.  Cloud  level  temperatures  at  various  grid  points 
during  SESAME  I. 
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Z.4.  Temperature  thresholds  of  satellite  imagery  for 
MLHR3 . 
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and  20 *  middle  cloud,  the  square  was  assumed  bo  be  covered 
with  70*  high  cloud.  It  was  then  assumed  that  latent  heat 
release  took  place  at  all  levels  below  the  level  of  cloud 
designated  as  the  predominant  cloud  in  a  grid  square.  If 
the  different  cloud  levels  covered  approximent 1 y  the  same 
area  of  a  grid  square,  the  highest  level  was  assumed  to  be 
the  predominate  cloud  type.  This  cloud  type  was  then 
designated  as  the  type  for  the  total  percentage  of  grid 
space  covered  by  cloud.  For  example,  if  a  grid  square  was 
covered  with  30*  low  cloud  and  30*  middle  cloud,  the 
square  was  then  assumed  to  be  covered  with  60*  middle 
cloud.  With  this  value  set,  the  modification  of  latent 
heating  values  was  done  using  the  same  procedure  as 
described  in  the  Lin  and  Smith  section.  As  before,  a  zero 
latent  heat  release  value  was  assigned  to  squares  with 
less  than  20*  of  their  areas  covered  with  cloud.  Also,  if 
clouds  were  observed  but  no  latent  heat  was  calculated, 
the  inverse  distance  weighted  interpolation  was  done  at 
each  level  using  the  initial  latent  heat  release  values. 
The  predominant  cloud  type  was  then  designated  and  the 
modification  was  applied.  Note  that,  as  with  the  original 
Lin  and  Smith  procedure(MLHRl ) ,  this  modification  was 
applied  to  the  TLHR ,  not  to  the  two  subcomponents 
seperately.  The  latent  heat  release  values  resulting  from 
this  modification  will  be  referred  to  in  the  results 
chapters  as  modified  latent  heat  release  thr ee ( MLHR3 ) . 


2.5  Evaluation  orocedures 


Initially,  latent  heat  release  values  were  evaluated  by 
comparing  the  areal  coverage  of  calculated  latent  heat 
release  with  radar  summaries  obtained  from  the  National 
Weather  Service,  satellite  cloud  depictions,  and  observed 
precipitation  fields.  Additional  comparisons  were  made  by 
calculating  area  averages  and  correlation  coefficients 
between  the  observed  precipitation  rates,  after  space  and 
time  smoothing,  and  the  precipitation  rates  corresponding 
to  the  vertically-integrated  values  of  calculated  total 
latent  heat  release.  Finally,  areal  latent  heating  and 
observed  precipitation  coverages  over  the  total  domain,  as 
well  as  homogeneous  sub-regions,  were  compared  by  simply 
counting  the  number  of  grid  points  over  which 
precipitation  was  observed  or  LHR  was  calculated.  All 
grid  point  observed  precipitation  and  calculated  LHR 
values  less  than  0.04mmh_1  were  regarded  as  trace  amounts 


and  were  set  to  zero. 


3. 


SYNOPTIC  DISCUSSION 


The  basic  synoptic  features,  with  corresponding 
precipitation  distributions,  are  represented  by  the  radar 
depictions  and  superimposed  surface  pressure  centers  and 
fronts  in  Fig. 3.1.  More  detailed  synoptic  discussions 
have  been  presented  for  the  SESAME  I  period  by  Alberty  et 
a  1 .  (1979),  Carlson  et  a  1  ■  (1980),  Moore  and  Fuelberg 
(1981),  and  Vincent  and  Homan  (1983).  The  following 
serves  as  a  brief  summary  of  the  evolution  of  the  cyclone 
and  its  associated  precipitation  field. 

Synoptic  conditions  at  1200GMT  10  April  featured  a  deep 
cyclone  located  over  northcentral  Colorado,  with  a  cold 
front  stretching  southward  through  New  Mexico.  During  the 
10th,  moist  air  from  the  Gulf  of  Mexico  was  flowing 
northward  across  Texas,  while  dry  air  from  eastern  New 
Mexico  pushed  a  dry  line  east  and  northeast  across  western 
Texas.  Early  radar  summaries  show  that  rainshower  and 
thundershower  activity  was  confined  to  the  region  affected 
by  a  weakening  upper  air  shortwave  that  was  moving  through 
the  Mississippi  Valley.  As  the  dry  line  moved  into 


westcentral  Texas  near  1800GMT,  there  was  a  dramatic 
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4.  THE  BASIC  LATENT  HEAT  RELEASE  CALCULATIONS 


V* 


4 . 1  Horizontal  distributions 

As  an  initial  step  in  the  evaluation  of  the  latent  heat 
calculations,  the  horizontal  fields  of  basic  latent  heat 
re  1  ease (BLHR)  estimates,  derived  from  the  parameterization 
schemes  described  in  Chapter  2,  are  examined.  Observed 
precipitation,  sate  1 1 i te— der i ved  cloud,  radar,  and  BLHR 
fields  for  four  map  times  are  presented  in  Figs . 4 . 1-4 . 4 . 
Before  proceeding  with  the  discussion,  it  must  be 
acknowledged  that  comparisons  of  these  data  forms  are 
complicated  by  a  measure  of  incompatabi 1 ity  that  exists 
between  them.  The  radar  maps  depict  instantaneous  fields 
of  larger  water  droplets,  some  of  which  are  associated 
with  subgrid-scale  events.  Satellite  imagery  allows  an 
instantaneous  view  of  cloud  tops,  but  obscures  the 
condensate  that  might  or  might  not  exist  below.  As  with 
the  former  two,  BLHR  values  are  obtained  at  particular 
times  but  are  representative  of  grid-scale  processes. 
Finally,  the  observed  precipitation  values  are  smoothed 
over  three-hour  periods(see  Sect. 2. 2)  to  make  them  more 
compatible  with  the  BLHR  values  and  are  subject  to  errors 
introduced  by  the  use  of  different  types  of  raingauges  and 
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different  precipitation  reporting  procedures.  Additional 
comment  on  raingauge  density  and  location  and  its  effect 
on  area  precipitation  totals  can  be  found  in  Woodley  £t 
al.  (1975)  and  Hildebrand  et  al .  (1979). 


Inconsistencies  between  the  radar  and  observed 
precipitation  summaries  are  particularly  troublesome 
because  each  represents  an  independent  standard  against 
which  the  computed  latent  heat  release  values  can  be 
compared.  Several  examples  of  such  inconsistencies  can  be 
identified  for  this  case  study.  At  1200GMT  (Fig.4.1),  an 
area  of  radar  echoes  in  northeastern  Texas  does  not  appear 
on  the  observed  precipitation  summary,  while  observed 
precipitation  areas  in  central  Texas  are  not  reflected  in 
the  radar  summary.  At  1800GMT  (Fig. 4. 2),  an  observed 
precipitation  pattern  extending  from  central  Texas  well 
into  eastern  Oklahoma  and  western  Arkansas  is 
substantially  larger  than  the  radar  area  which  is  confined 
to  northern  Texas.  In  this  latter  example  the 
precipitation  area  occurs  in  response  to  a  rapidly 
expanding  area  of  convection  which  spread  north  and  east 
out  of  northcentral  Texas.  Interestingly,  while  the 
instantaneous  radar  summary  does  not  depict  this  area,  a 
combination  of  successive  radar  summaries  does.  This  is 
clearly  shown  in  Fig. 4. 5,  which  depicts  the  radar 
summaries  for  1735GMT,  1835GMT,  and  1935GMT.  This  figure 
also  shows  yet  another  inconsistency  example,  namely  radar 
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echoes  in  eastern  Texas  and  western  Louisiana  which  are 
not  reflected  in  the  observed  precipitation  summary. 
Unfortunately,  the  failures  of  the  radar  and  observed 
precipitation  summaries  to  match  do  not  follow  a 
systematic  pattern.  Therefore,  it  will  be  necessary  to 
note  the  specific  impacts  of  these  data  differences  as 
they  become  apparent  in  the  qualitiative  comparisons  of 
the  LHR  fields.  Despite  these  differences,  it  is  hoped 
that  taken  together  these  data  forms  can  yield  useful 
insight  into  the  quality  of  the  BLHR  calculations,  as  well 
as  the  impact  of  the  three  modification  schemes. 


4.1.1  1200GMT  April  10,1979 

The  BLHR  field  at  this  time  (Fig. 4.1)  compares 
favorably  with  the  other  data  forms  in  several  areas.  The 
arc  of  observed  precipitation  which  appears  as  a  broken 
pattern  curving  from  eastern  Colorado  through  Kansas  and 
Missouri  and  into  Arkansas  is  generally  matched  by  both 
the  satellite-derived  cloud  cover  and  the  radar  summary. 
BLHR  values  also  reflect  this  precipitation  pattern  even 
to  the  extent  of  pinpointing  the  observed  precipitation 
maximum  in  eastern  Arkansas.  In  addition,  the  BLHR 
maximum  in  eastern  Texas  can  be  justified  by  the  radar 
reports  near  the  Louisiana  border,  although  this  rapidly 
moving  area  of  convection  is  not  reflected  in  the  observed 
precipitation  field.  Also,  observed  precipitation  in 


40 


northcentral  Texas  extending  into  Oklahoma  is  matched  by 
the  calculated  precipitation  values.  A  noteable  failure 
of  the  BLHR  calculations  is  apparent  in  Mississippi  where 
the  presence  of  relatively  strong  BLHR  values  is  not 
supported  by  either  the  radar  or  the  observed 
precipitation  summaries. 


£* 


4.1.2  1800GMT  April  10,1979 

The  BLHR  calculations  at  this  time  (Fig. 4. 2)  indicate  a 
large  area  of  latent  heat  release  extending  from  New 
Mexico  northeastward  through  the  panhandle  of  Texas  into 
Kansas  and  Missouri.  Satellite  imagery  verifys  the 
presence  of  clouds  in  this  region,  and  the  radar  summary 
shows  good  agreement.  However,  the  radar  echoes  are 
scattered  and  the  tops  are  not  very  high,  features  that 
perhaps  explain  the  lack  of  observed  precipitation  in 
portions  of  this  region.  BLHR  calculations  perform  well 
in  northcentral  Texas  as  they  correctly  identify  an  area 
of  deep  convection  which  is  shown  on  the  radar  summary. 
The  observed  precipitation  field  is  also  consistent  with 
these  results.  As  previously  noted  (Fig. 4. 5),  this  area 
of  convection  grew  rapidly  in  the  next  two  hours  expanding 
in  all  directions.  This  perhaps  explains  why  the  observed 
precipitation  area  extends  well  into  eastern  Oklahoma  and 
western  Arkansas.  The  BLHR  summary  also  indicates  an  area 


of  latent  heating  in  southeastern  Texas  and  Louisiana. 


The  radar  summary  is  in  good  agreement,  and  the  satellite 
depiction  also  shows  extensive  cloud  cover  in  this  area. 
However,  this  region  failed  to  record  any  precipitation. 

The  BLHR  values  in  northern  Missouri  are  not  reflected 
in  the  radar  fields,  nor  by  the  observed  precipitation 
summary.  However,  cloud  cover  is  indicated  in  the 
satellite  depiction.  In  this  case  the  BLHR  calculations 
were  dominated  by  the  SLHR  component.  Further,  stratiform 
clouds  are  indicated  in  the  satellite  photographs.  Thus, 
this  may  be  an  area  of  non-precipitating  stratiform  clouds 
yielding  echoes  too  weak  for  radar  identification. 

4.1.3  OOOOGMT  April  11,1979 

The  BLHR  estimates  for  this  time  (Fig. 4. 3)  show  a  large 
area  of  latent  heating  throughout  much  of  the  SESAME  I 
grid,  with  the  primary  latent  heating  axis  extending  from 
southern  Missouri  through  central  Oklahoma  and  into 
northcentral  Texas.  The  radar  and  observed  precipitation 
summaries  concur  with  this.  The  BLHR  calculations  also  do 
well  in  pinpointing  the  observed  precipitation  noted  in 
Louisiana  by  the  other  summaries.  In  addition,  the 
absence  of  observed  precipitation  in  Arkansas  and 
southwestern  Texas  is  similarly  depicted  in  the  BLHR 
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The  BLHR  estimates  are  less  successful  in  northeastern 
Texas  and  southeastern  Oklahoma,  where  calculations 
suggest  latent  heating  which  is  not  confirmed  by  any  of 
the  other  three  data  forms.  The  BLHR  summary  also  misses 
the  important  region  of  convection  south  of  the  Texas- 
Oklahoma  border  near  the  southwestern  corner  of  Oklahoma. 
All  the  other  summaries  indicate  its  presence. 


•  V 


4.1.4  0600GMT  April  11,1979 

At  this  final  map  time  (Fig. 4. 4),  the  primary  BLHR 
region  has  moved  eastward  with  its  orientation  continuing 
along  a  northeast-southwest  axis.  Both  the  radar  and  the 
observed  precipitation  summaries  agree  that  the  majority 
of  the  precipitation  extends  from  northeast  Missouri 
through  eastern  Oklahoma  to  central  Texas.  The  BLHR  field 
captures  the  northern  two-thirds  of  this  area.  The 
observed  precipitaion  maxima  are  mirrored  by  BLHR  maxima 
in  northeastern  Missouri  and  northeastern  Oklahoma.  Also, 
a  secondary  area  of  convection  in  western  Kansas  shows  up 
well  on  all  the  summaries,  as  does  a  region  of  snow 
showers  in  eastern  Colorado. 

The  BLHR  field  fails  at  two  noteable  locations.  The 
first  is  the  severe  convection  seen  in  the  radar  summary 
in  northcentral  Texas.  The  BLHR  calculations  for  this 
area  show  no  convective  component  of  latent  heat  release 
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and  vary  small  SLHR  values.  The  second  is  the  BLHR 
precipitation  field  in  eastern  Texas  and  western 
Louisiana.  The  radar  summary  shows  that  an  important 
radar  report  in  southwestern  Louisiana  is  missing.  Still, 
both  the  satellite-derived  cloud  cover  and  the  observed 
precipitation  summaries  show  no  clouds  and  no 
precipitation  reports  for  this  area. 


4.1.5  Synopsis  of  horizontal  depictions 

The  overview  presented  here  shows  that  the  BLHR  fields 
are  reasonably  successful  in  identifying  the  major  areas 


of  precipitation. 

part i eu 1 ar 1 y 

after 

the 

convection 

becomes  organized. 

The  failure  of 

the 

parameterzation 

scheme  to  compare  favorably 

with 

the 

observed 

precipitation  fields  in  some  areas  can  be  attributed  to  a 
variety  of  factors.  The  first  derives  from  the  fact  that 
parameterization  schemes  are  intended  to  depict  grid-scale 
convective  and  stable  latent  heating  fields.  In  fact,  the 
observed  precipitation  and  radar  summaries  show  the 
precipitation  to  be  largely  subgrid-scale  through  1800GMT, 
10  April.  By  2100GMT  the  synoptic-scale  forcing  becomes 
more  dominant,  yielding  more  organized  precipitation 
distributions.  As  a  result,  the  BLHR  values  are  more 
comparable  with  observed  precipitation  values  at  these 
later  times.  The  contrast  between  the  comparability  of 


these  two  fields  at  the  earlier  and  later  times  will  be 


and  radar  summaries;  thus,  the  latter  are  of  little  value 


m 
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as  comparative  parameters.  Instead,  the  vertical 
distributions  of  latent  heat  release  for  the  present  study 
are  compared  to  those  from  two  other  studies,  the  vertical 
profiles  of  area-averaged  heating  rates  by  Fuelberg  et 
a  1 . ( 1 985 )  and  time/height  cross-sections  by 
Robertson ( 1983) .  In  all  cases  the  quantities  are  area 
averaged  over  the  entire  computational  domain. 


4.2.1  Vertical  distributions  of  CLHR 

The  rapid  increase  in  convective  activity  during  the 
period  under  consideration  is  readily  seen  in  the  radar 
summaries  presented  in  Figs . 4 . 1-4 . 4 .  This  increased 
convection  is  also  indicated  by  the  marked  increase  in 
CLHR  that  took  place  after  1500GMT  10  April  (Fig. 4. 6). 
The  maximum  CLHR  occurred  in  the  400  mb-500  mb  layer  at 
0000GMT. 


When  these  cross-sections  are  compared  with  the 
vertical  profiles  of  Fuelberg  et  a  I .  ( 1985) (Fig . 4 . 7) , 

|  generally  good  agreement  is  seen.  Both  distributions 

I 

indicate  increasing  CLHR  through  0000GMT  11  April  and  then 
decreasing  at  0600GMT,  with  the  weakest  CLHR  occurring  at 
j  1200GMT  and  the  strongest  at  0000GMT.  This  was  true  at 

i 

nearly  every  pressure  level.  Also,  Fuelberg's  vertical 
profiles  agree  that  the  maximum  CLHR  values  occurred  in 
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Figure  4.6.  Time/height  cross-section  of  BLHR  convective 
component(deg  day'1). 
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the  400-500  mb  layer  for  the  latter  three  time  periods  and 
that  the  bulk  of  the  CLHR  took  place  between  800  mb  and 
200  mb.  Similarities  in  the  two  sets  of  results  are  not 
surprising  since  Fuelberg's  computational  schemes  were 
similar  to  those  used  in  this  study.  Differences  that  are 
apparent  in  the  two  sets  of  results  can  be  attributed  to 
Fuelberg’s  use  of  a  different  data  analysis  procedure,  a 
somewhat  larger  grid  spacing,  and  slightly  different 
modifications  of  the  basic  parameterization  scheme.  In 
particular,  it  should  be  noted  that  most  of  Fuelberg’s 
values  ranged  from  2-80Cday_1,  while  the  results  of  this 
study  were  a  factor  of  two  higher. 

A  CLHR  vertical  cross  section  was  also  obtained  from 
the  work  done  by  Robertson (1983) (Fig . 4 . 8) .  In  his  study, 
the  grid-scale  heating  from  convection  was  derived  using  a 
mass  flux  technique  which  incorporated  a  1-D  cloud  model 
and  includes  consideration  of  precipitation  observed  at 
the  ground  and  the  cloud  structure  observed  from  SMS-I/ 
GOES  IR  imagery.  Since  Robertson’s  scheme  does  not  depend 
on  the  Kuo  parameterization,  it  yields  much  more 
independent  comparative  statistics  than  does  the  Fuelberg 
et  a  1 .  (1985)  paper.  In  general,  Robertson's  results 
contain  the  same  character istics  found  in  the  vertical 
CLHR  distributions  of  this  study.  However,  Robertson's 
area-averaged  CLHR  values  were,  on  the  whole,  a  factor  of 
two  smaller  than  those  calculated  in  this  study. 


ROBERTSON'S  CONV.  HEATING  RATE 
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4.2.2  Vertical  distributions  of  SLKR 

The  SLHR  vertical  distributions(Fig . 4. 9)  also  show  an 
increase  throughout  the  period  of  this  study.  This  is  in 
response  to  an  increased  flow  of  moist  air  from  the  Gulf 
of  Mexico  over  the  eastern  two-thirds  of  the  area.  This, 
coupled  with  an  expanding  area  of  upward  motion  (Moore  and 
Fuelberg,  1981).  supports  the  larger  SLHR  values  found 
late  in  the  period.  When  compared  with  the  Fuelberg  et 
a  1 .  (1985)  vertical  prof i 1 es (Fig . 4 . 1 0) ,  both  show  a 
general  increase  of  SLHR  values  at  nearly  all  levels 
throughout  the  period.  However,  as  with  the  CLHR  results, 
the  SLHR  values  in  this  study  are  on  the  whole  larger  than 
those  found  in  Fuelberg.  The  variations  mentioned  earlier 
(see  section  4.2),  as  well  as  his  use  of  90%  relative 
humidity  as  a  discriminator  for  SLHR  calculations 
(compared  to  the  80%  used  in  this  study),  can  be 
identified  as  reasons  for  differences  seen  in  the  two  sets 
of  results. 
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5.  THE  MODIFIED  LATENT  HEAT  RELEASE  CALCULATIONS 

5 . 1  Horizontal  distributions 

This  chapter  presents  an  analysis  of  the  changes  in  the 
latent  heat  release  fields  brought  about  by  the  three 
modification  schemes  summarized  in  chapter  2.  As  in  the 
preceding  chapter  comparisons  are  made  between  the 
modified  latent  heat  release  f ie 1 ds (MLHR1 , 2 , 3) , radar 
patterns,  and  observed  precipitation  summaries.  These 
fields  and  the  previously  discussed  BLHR  values  are 
presented  in  Figs . 5. 1-5. 4. 


5.1.1  1200GMT  April  10,  1979 

The  unenhanced  satellite  cloud  image (Fig . 4 . 1 )  for  this 
time  results  in  the  retention  of  the  arc  of  BLHR 
values (Fig . 5 . 1 )  that  extends  from  eastern  Colorado 
eastward  through  Kansas  and  Missouri  and  then  southward 
through  Arkansas.  This  area,  which  is  also  reflected  in 
the  observed  precipitation  and  radar  summaries,  is 
presented  as  nearly  continuous  in  the  MLHR1  and  MLHR3 
fields.  The  MLHR2  modification  presents  a  much  more 


broken  pattern  depiction,  which  matches  more  closely  the 


Figure  5. 1 


Top  left:  NMC  radar  summary  chart  for  1200GMT 
10  April  1979.  Cloud  tops  in  dam,  motion 
vectors  in  m  s"1. 

Top  right:  Observed  precipitation  rate(mm  h_1) 
for  1200GMT  10  April  1979. 

Middle  left:  Vertically  integrated  BLHR 
combining  convective  and  stable  components (mm 
h*1). 

Middle  right:  Vertically  integrated  MLHR1 
combining  convective  and  stable  components (mm 
h->). 

Bottom  left:  Vertically  integrated  MLHR2 
combining  convective  and  stable  components(mm 
h->>. 

Bottom  right:  Vertically  integrated  MLHR3 
combining  convective  and  stable  components (mm 
h 
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Figure  5.2.  Top  left:  NMC  radar  summary  chart  for  1800GMT 
10  April  1979.  Cloud  tops  in  dam,  motion 
vectors  in  m  s"1. 

Top  right:  Observed  precipitation  rate(mm  h"1) 
for  1800GMT  10  April  1979. 

Middle  left:  Vertically  integrated  BLHR 
combining  convective  and  stable  components (mm 
h-1). 

Middle  right:  Vertically  integrated  MLHR1 
combining  convective  and  stable  components (mm 
h-‘>. 

Bottom  left:  Vertically  integrated  MLHR2 
combining  convective  and  stable  components(mm 
h-«). 

Bottom  right:  Vertically  integrated  MLHR3 
combining  convective  and  stable  components(mm 
h  ~  1  )  . 


Figure  5.3.  Top  left:  NMC  radar  summary  chart  for  OOOOGMT 
11  April  1979.  Cloud  tops  in  dam,  motion 
vectors  in  m  s" 1 . 

'  Top  right:  Observed  precipitation  rate(mm  h"1) 
for  OOOOGMT  11  April  1979. 

Middle  left:  Vertically  integrated  BLHR 
combining  convective  and  stable  components (mm 
h-1). 

Middle  right:  Vertically  integrated  MLHR1 
combining  convective  and  stable  components (mm 
h-‘). 

Bottom  left:  Vertically  integrated  MLHR2 
combining  convective  and  stable  components (mm 
h-1). 

Bottom  right:  Vertically  integrated  MLHR3 
combining  convective  and  stable  components ( (mm 
h-O. 
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Figure  5.4. 


Top  left:  NMC  radar  summary  chart  for  0600GMT 
11  April  1979.  Cloud  tops  in  dam,  motion 
vectors  in  m  s “ 1 . 

Top  right:  Observed  precipitation  rate(mm  h" 1 ] 
for  0800GMT  11  April  1979. 

Middle  left:  Vertically  integrated  BLHR 
combining  convective  and  stable  components (mm 
h->). 

Middle  right:  Vertically  integrated  MLHR1 
combining  convective  and  stable  components (mm 
h->). 

Bottom  left:  Vertically  integrated  MLHR2 
combining  convective  and  stable  components (mm 
h-‘). 

Bottom  right:  Vertically  integrated  MLHR3 
combining  convective  and  stable  components (mm 
h-»>. 
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general  characteristics  of  the  observed  precipitation 
pattern.  A  similar  result  is  seen  over  New  Mexico,  where 
the  MLHR1  and  MLHR3  schemes  increase  the  areal  coverage  of 
calculated  latent  heat  release.  The  MLHR2  pattern  again 
suggests  spotty  latent  heating,  in  general  agreement  with 
both  the  observed  precipitation  and  radar  fields. 

5.1.2  1800GMT  April  10.  1979 

The  principle  region  of  latent  heat  release  extending 
from  New  Mexico  northward  through  Colorado  and  southern 
Nebraska  is  preserved  by  all  three  modifications,  although 
it  is  perhaps  too  small  in  the  MLHR2  depiction,  Further,  a 
small  latent  heat  release  area  in  southern  Louisiana  is 
identified  in  all  of  the  modifications.  This  feature  is 
supported  by  the  radar  summary  but  is  not  depicted  in  the 
observed  precipitation  summary. 

The  most  striking  differences  are  seen  in  northern 
Texas,  Oklahoma,  Arkansas,  and  Missouri.  Both  the  MLHR1 
and  MLHR3  modifications  show  extensive  latent  heating  in 
these  areas,  much  of  which  is  unsupported  either  by  the 
radar  or  the  observed  precipitation  summaries.  The  MLHR2 
field  presents  a  much  more  satisfactory  pattern,  although 
it  clearly  underestimates  the  extent  of  the  observed 
precipitation  in  northern  Texas  and  Oklahoma.  All  of  the 
modifications  fail  to  extend  latent  heating  values  into 
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central  Texas  to  account  for  the  precipitation  that  is 
indicated  by  the  observed  precipitation  summary. 

5.1.3  OOOOGMT  April  11.  1979 

An  extensive  precipitation  pattern  depicted  by  both  the 
radar  and  the  observed  precipitation  summaries  extends 
from  northern  Texas  through  Oklahoma  into  Kansas,  and  then 
eastward  into  Missouri  and  Arkansas.  All  modifications 
reflect  this  region  very  well,  even  to  the  extent  of 
identifying  the  axis  of  heaviest  precipitation  which 
extends  from  southern  Missouri  through  central  Oklahoma 
and  into  northcentral  Texas.  Both  the  MLHR1  and  MLHR3 
modifications  indicate  latent  heat  release  in  Louisiana, 
separating  this  area  from  the  larger  precipitation  region 
to  the  north.  The  MLHR2  pattern  also  eliminates  latent 
heating  in  eastern  Texas  and  southwestern  Arkansas,  but 
its  Louisiana  values  are  not  as  well  placed  as  the  other 
modifications. 


5.1.4  0600GMT  April  11,  1979 

Once  again,  all  of  the  modifications  mirror  the  large 
precipitation  area,  including  the  maximum  precipitation 
axis  extending  from  northeastern  Missouri  through  eastern 
Oklahoma  and  into  central  Texas.  The  modifications  are 


able  to  pinpoint  the  northern  two  precipitation  maxima, 
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one  in  northeastern  Missouri  and  the  other  in  northeastern 
Oklahoma,  and  also  reflect  the  relatively  unbroken  nature 
of  this  axis  as  it  extends  into  central  Texas.  This 
pattern  appears  to  be  more  reasonable  than  the  broken 
pattern  of  latent  heat  release  suggested  by  the  BLHR 
field.  In  addition,  the  MLHR1  and  MLHR3  modifications  in 
this  area  extend  values  further  south  into  Texas  than  is 
seen  in  the  BLHR  and  MLHR2  calculations,  and  all  three 
modifications  successfully  eliminate  BLHR  values 
calculated  in  eastern  Texas  and  western  Louisiana. 
Finally,  contrary  to  the  other  two  modifications,  the 
MLHR3  field  extends  into  the  Texas  panhandle  and  northern 
New  Mexico,  a  feature  that  is  not  reflected  by  the  radar 
and  observed  precipitation  summaries. 


5.1.5  Synopsis  of  horizontal  depictions 


In  evaluating  the  various  modifications  made  to  the 
BLHR  calculations,  some  important  points  become  readily 
apparent.  In  cases  where  the  BLHR  calculations  did  very 
poorly  in  reflecting  observed  precipitation,  the 
modifications  were  able  to  make  some  dramatic  changes  in 
the  areal  coverage  of  the  calculated  precipitation  fields. 
This  was  found  to  be  true  especially  early  in  the  study 
period.  After  2100GMT,  the  BLHR  calculations  were  much 
more  successful  in  identifying  the  pattern  of  observed 
precipitation.  At  these  later  times  the  modifications 
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made  changes  in  some  areas,  but  the  overall  latent  heating 
pattern  remained  relatively  unaffected.  It  can  also  be 
seen  that  neither  the  BLHR  calculations  nor  the 
modifications  succeeded  in  capturing  isolated  convective 
activity  well.  This  isolated  convection  was  most 
prevalent  in  the  two  early  time  periods  of  1200GMT  and 
1800GMT . 

It  is  premature  at  this  point  to  advance  one 
modification  as  superior  to  the  others  or  to  state 
conclusively  that  any  of  them  are  more  precise  than  the 
BLHR  calculations.  Yet  several  important  observations  can 
be  made  concerning  the  different  modifications  used  in 
this  study.  The  calculated  precipitation  fields  produced 
by  the  235K  enhancement  modification  (MLHR2)  consistently 
matched  the  radar  summary  much  closer  than  the  fields  of 
the  other  modifications.  This  supports  the  contention 
made  in  a  study  done  by  Arkin(1979)  that  the  235K 
threshold  is  able  to  isolate  the  regions  of  convective 
precipitation  occurring  in  a  cloud  field.  However,  the 
MLHR2  modification  also  failed  to  identify  some  areas  of 
light  precipitation  which  were  captured  by  the  raingauge 
network.  While  the  235K  enhancement  reduced  the  latent 
heating  estimate  by  decreasing  the  areal  coverage  of 
latent  heat  release,  the  3-layer  enhancement (MLHR3) 


typically  increased  the  areal  coverage  of  calculated 
latent  heating.  This  was  accomplished  by  producing  values 
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at  points  that  were  missed  in  the  BLHR  calculations,  but 
where  there  was  cloud  cover.  While  under  certain 
circumstances  this  caused  the  calculated  precipitation 
patterns  to  better  reflect  the  radar  and  observed 
precipitation  summaries,  improvement  was  not  always  noted. 
Frequently,  this  third  modification  returned  latent  heat 
release  values  to  regions  which  had  been  vacated  by  the 
main  precipitation  area,  but  in  which  cloud  cover  lingered 
behind.  Finally,  the  MLHR1  and  MLHR3  modifications 
calculated  latent  heat  release  in  regions  that  were 
influenced  by  the  cloud  cover  preceding  and  surrounding  a 
major  precipitation  region.  This  tendency  for  the 
calculated  latent  heat  release  fields  to  envelope  the 
observed  precipitation  fields  matches  some  of  the 
conclusions  found  in  the  work  done  by  Lin  and  Smith(1985) 
and  may  reflect  the  ability  of  the  computational  schemes 
to  identify  latent  heat  release  in  non-precipitating 
regions . 


5 . 2  Vertical  cross-sections 

The  modifications  used  in  this  study  affect  the 
vertical  distribution  of  th<  BLHR  values  in  different 
ways.  As  described  in  Chapter  2,  the  MLHR1  and  MLHR3 
modifications  impact  the  total  BLHR  values,  while  the 
MLHR2  scheme  modifies  the  two  components  of  latent  heat 
release  separately.  Therefore,  a  complete  evaluation  of 
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the  vertical  latent  heat  release  distributions  is  limited 
to  the  TLHR  values  for  both  the  basic  and  modified 
calculations.  These  in  turn  can  be  compared  with  vertical 
profiles  of  the  sum  of  the  two  components  of  LHR  presented 
by  Fuelberg  et  a  1 .  (1985)  (see  section  4.2). 

Vertical  cross-sections  of  BLHR,  MLHR1 ,  MLHR2 ,  and 
MLHR3  are  presented  in  Fig. 5. 5.  Several  important  points 
are  immediately  apparent.  This  figure  shows  that  the 
basic  pattern  of  BLHR  values  remains  relatively  unaltered 
by  the  modifications.  All  cross-sections  show  latent  heat 
release  increasing  at  virtually  all  levels  prior  to 
0300GMT,  with  the  bulk  of  LHR  occurring  between  750  mb  and 
350  mb.  Also,  each  modification  reflects  the  double 
maxima  seen  near  700  mb  and  400  mb  at  0300GMT  in  the  BLHR 
calculation. 

The  greatest  impact  from  the  modifications  is  on  the 
magnitudes  of  the  latent  heat  release  values.  The  MLHR2, 
and  to  a  lesser  extent  the  MLHR1  scheme,  reduces  the  area- 
averaged  total  latent  heat  release  values  at  nearly  every 
time  and  level,  while  the  MLHR3  modification  retains 
nearly  the  same  magnitudes.  These  conclusions  are  similar 
to  the  observations  presented  in  the  evaluation  of  the 
horizontal  LHR  fields;  that  is,  the  MLHR2  modification  has 
the  largest  impact  on  the  BLHR  values,  and,  in  fact, 
reduces  the  latent  heating  estimates.  Also,  while  the 
MLHR3  modification  succeeds  in  eliminating  LHR  values  at 


Figure  5.5.  Top  left:  Vertical  cross-section  of  BLHR 

combining  convective  and  stable  components(deg 
day” 1 ) . 

Top  right:  Vertical  cross-section  of  MLHR1 
combining  convective  and  stable  components (deg 
day" 1 ) . 

Bottom  left:  Vertical  cross-section  of  MLHR2 
combining  convective  and  stable  components(deg 
day- 1 ) . 

Bottom  right:  Vertical  cross-section  of  MLHR3 
combining  convective  and  stable  components(deg 
day" 1 ) . 


some  gr*i<l  points  and  at  various  levels  in  the  vertical, 
this  modification  also  returns  values  at  other  points  and 
levels,  thus  leaving  the  overall  average  values  relatively 
unaffected . 


The  vertical  profiles  of  the  two  components  of  LHR 
provided  in  the  study  of  Fuelberg  et  a  1 .  (1985)  were 
summed  at  each  pressure  level  in  order  to  obtain  a 
vertical  distribution  of  TLHR.  The  resulting  profile  is 
shown  in  Fig. 5. 6.  Fuelberg  provided  values  for  only  four 
time  periods  that  were  common  with  the  times  used  in  this 
study.  Fuelberg’s  calculations  indicate  that  the  maximum 
latent  heating  takes  place  between  OOOOGMT  and  0600GMT,  in 
good  agreement  with  this  study's  results.  Both 
distributions  have  the  bulk  of  the  LHR  occurring  between 
700  mb  and  350  mb.  For  the  latter  three  time  periods 
Fuelberg  indicated  that  the  maximum  LHR  occurred  between 
400  mb  and  550  mb,  also  in  agreement  with  this  study’s 
findings.  Of  the  three  modifications,  the  MLHR2  vertical 
cross-section  provides  LHR  values  that  are  closer  in 
magnitude  to  Fuelberg ’s  results.  Differences  in  the 
vertical  distributions  of  LHR  found  in  the  modified 
calculations  and  those  of  Fuelberg  are  similar  to  those 
found  in  the  BLHR  estimates. 


Finally,  since  the  MLHR2  results  are  the  most 
comparable  with  the  Fuelberg  et  a  I .  (1985)  profiles,  it  is 
of  interest  to  compare  the  convective  component  of  these 


results  with  those  of  Robertson ( 1983) ,  discussed  in 
section  4.2.  Fig. 5. 7  shows  the  cross-sections  of  the 
MLHR2  convective  component  together  with  that  of 
Robertson.  Comparing  this  with  the  BLHR  convective 
component  displayed  in  Fig. 4. 6,  it  is  clear  that,  although 
still  larger,  the  MLHR2  values  are  significantly  closer  to 
the  Robertson's  results  than  are  the  basic  convective 


Ml_HR  2  iso  ROBERTSON'S  CONV.  HEATING  RATE 


of  convective 


6.  FURTHER  EVALUATION  OF  BASIC  AND  MODIFIED  LHR 

CALCULATIONS 


6 . 1  Area-averaged  values 

Despite  the  shortcomings  mentioned  in  section  4.1 
concerning  the  use  of  gridded  precipitation  totals,  these 
values  remain  the  most  independent  standard  available  for 
an  evaluation  of  LHR  calculations.  This  evaluation  serves 
to  compliment  the  conclusions  drawn  earlier  in  the 
qualitative  comparison  of  horizontal  and  vertical  LHR 
distributions.  The  simplest  and  least  detailed  statistics 
which  can  be  considered  in  this  evaluation  are  the  area- 
averaged  observed  precipitation  and  the  basic  and  modified 
latent  heating  values,  summarized  in  Fig. 6.1.  This  figure 
supports  the  earlier  qualitative  observation  that 
precipitation  rates  increased  over  the  SESAME  I  grid  from 
1200GMT  10  April  to  0300GMT  11  April.  The  rapid  growth  of 
convection  that  was  depicted  in  the  radar  summaries  of 
Fig. 4. 5  is  reflected  by  the  large  increase  in  all  of  the 
area-averaged  quantities  after  1800GMT.  This  period  of 
large  increase  corresponds  to  the  time  interval  during 


which  synopt i c— sea  1 e  forcing  became  more  dominant  and  more 
organized  precipitation  distributions  were  evident(see 


Section  4.1.5).  The  rapid  increase  in  precipitation 
amount  and  its  relationship  to  synoptic— sea  1 e  forcing 
suggest  that  it  would  be  appropriate  to  summarize  these 
and  later  statistics  by  averaging  over  two  periods. 
Period  1  encompasses  the  first  three  map  times  (1200GMT, 
1500GMT,  1800GMT),  while  Period  2  encompasses  the  last 
four  map  times (2100GMT ,  OOOOGMT,  0300GMT,  0600GMT) . 


Fig. 6.1  shows  that,  with  the  exception  of  0600GMT,  all 
three  modifications  reduce  the  overestimate  of  latent 
heating  seen  in  the  BLHR  totals.  The  MLHR1  and  MLHR3 
calculations  remain  larger,  while  the  MLHR2  are  smaller, 
than  the  observed  precipitation  rates.  The  extent  of 
these  modifications  is  summarized  by  Table  6.1,  which 
presents  the  average  absolute  departures  from  the  area- 
averaged  observed  precipitation  totals  for  the  basic  and 
modified  LHR  values.  Overall,  the  MLHR1  modification 
yields  the  smallest  deviation  from  the  observed  area 
averaged  totals;  however,  in  the  first  period  the  MLHR2 
departure  is  the  smallest.  An  analysis  of  variance  using 
an  F— test  (Hicks, 1982)  was  done  on  the  mean  departures  of 
area-averaged  basic  and  modified  LHR  calculations  from  the 
area-averaged  observed  precipitation  values.  Using  the 
assumption  that  an  underestimation  of  area  averaged 
observed  precipitation  does  not  hold  any  advantage  over  an 
overest imat i on ,  there  was  found  to  be  no  statistical 


difference  between  the  mean  LHR  departures,  at  the  5% 
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Table  6.1.  Departure  of  area— averaged  LHR  quantities  from 
area— averaged  observed  precipitation  values. 

DEPARTURE 

FROM  OBSERVED 

PARAMETER 

OVERALL 

PERIOD  1 

PERIOD  2 

BLHR 

.21 

.19 

.24 

MLHR1 

.12 

.09 

.15 

MLHR2 

-.15 

-.06 

-.21 

MLHR3 

.15 

.14 

.16 

PERIOD  1 
PERIOD  2 

1200GMT  -  1800GMT 
2100GMT  -  0600GMT 

8 


level  of  significance.  Therefore,  these  statistics  alone 
offer  little  to  suppo-t  the  superiority  of  one  set  of 
calculations  over  the  others. 

6 . 2  Correlation  coefficients 

To  compare  the  overall  grid  point  distribution  of 
observed  precipitation  with  those  of  the  basic  and 
modified  LHR  values,  these  correlation  coefficients  were 
computed  and  are  shown  in  Fig, 6. 2.  The  ability  of  the 
various  latent  heat  release  calculations  to  reflect  the 
observed  precipitation  distribution  varies  greatly  from 
Period  1  to  Period  2.  At  the  first  three  map  times,  when 
the  precipitation  field  was  largely  subgrid-scale, 
correlations  are  very  small.  After  1800GMT  the 
correlation  coefficients  rise  rapidly  for  all  of  the 
schemes.  Fig. 6. 2  also  shows  that  at  every  map  time  all  of 
the  modifications  have  higher  correlation  coefficients 
than  do  the  BLHF  calculations,  with  the  MLHR2  modification 
yielding  the  best  correlations. 

Table  8.2  summarizes  the  Fig. 6. 2  results  by  focusing  on 
the  percentage  improvement  in  the  correlation  coefficients 
produced  by  the  modifications.  This  table  again  shows 
that  all  three  modifications  improved  the  BLHR 
calculations.  The  relative  impact  of  these  modifications 

1  , 


is  especially  noteable  in  Period 


a  1  though 
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Table  6.2.  Percent  improvement  of  modified  over  basic  LHR 
correlation  coef f icients . 


PERCENT  IMPROVEMENT 
FROM  BLHR  VALUES 


PERIOD  2  2100GMT  -  0600GMT 


%  Improvement  =  MOD  VALUE  -  BLHR 

BLHR  X  100,0 


MODIFICATION 

OVERALL 

PERIOD  1 

PERIOD  2 

MLHRI 

23% 

73% 

10% 

MLHR2 

45% 

101% 

28% 

MLHR3 

23% 

60% 

14% 

PERIOD  1 

1200GMT  -  1800GMT 

8 


correlations  remain  very  poor.  The  MLHR2  modification 
yields  the  greatest  improvement  in  both  periods.  An 
analysis  of  variance  using  an  F-test  (Hicks, 1982)  was  done 
on  the  mean  correlation  coefficients  for  the  basic  and  the 
modified  calculations.  It  was  found  that  there  was  a 
significant  variation  in  these  values.  A  Newman-Kuels 
range  test  was  done  to  help  pinpoint  where  the  variation 
in  mean  values  occurred.  This  test  showed  that  the  MLHR2 
mean  correlation  coefficient  was  significantly  different 
than  the  BLHR  value,  but  that  there  was  no  significant 
difference  between  the  BLHR  mean  value  and  the  MLHRl  and 
MLHR3  mean  values. 

6 -  3  Grid-point  precipitation/LHR  comparison 

The  correlation  coefficient  results  reveal  a 
fundamental  problem  with  this  statistic.  When  the 
precipitation  is  subgrid-scale  in  character,  it  is  very 
difficult  for  any  computational  scheme  which  uses  grid- 
scale  data  to  specify  LHR  in  sufficient  detail  to  yield 
high  correlations.  Thus,  while  various  schemes  may  yield 
LHR  in  appropriate  regions,  they  likely  will  not  match  the 
exact  locations  of  the  observed  precipitation.  As  a 
consequence,  even  though  from  a  dynamical  standpoint  a 
reasonable  LHR  distribution  is  presented,  the  correlation 
coefficient  may  false'y  indicate  that  the  calculations  are 
of  little  value.  This  is  indeed  the  situation  at  the 
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first  three  map  times,  when  Fig. 6. 2  shows  very  poor 
correlations  but  Figs. 5.1  and  5.2  suggest  a  more  favorable 
comparison . 

An  alternate  way  of  comparing  LHR  distributions  is  to 
determine  for  a  given  area  the  total  number  of  points  at 
which  precipitation  was  observed  and  LHR  was  calculated  at 
each  map  time  of  the  study.  This  compilation  for  the 
total  grid  area  is  given  in  Table  6.3.  These  totals  serve 
to  reinforce  some  earlier  conclusions  concerning  the 
changes  in  precipitation  and  LHR  patterns  with  time.  The 
general  expansion  in  these  areas  prior  to  0300GMT  is 
clearly  shown  by  the  increase  in  the  number  of  LHR  grid 
points  within  each  field.  The  tendency  for  the  MLHRl  and 
MLHR3  modifications  to  increase  the  area  of  LHR(see 
section  5.1.5)  is  shown  by  the  larger  number  of  points 
where  latent  heating  is  calculated  when  compared  with  the 
BLHR  values.  The  decrease  in  the  calculated  latent  heat 
release  area  caused  by  the  MLHR2  scheme  is  also  seen  in 
Table  6.3.  Clearly,  gridpoint  totals  for  the  MLHR2  scheme 
compare  more  closely  with  the  observed  precipitation 
values  than  do  the  other  computational  schemes. 

From  a  dynamical  standpoint  it  is  sufficient  to 
properly  specify  the  LHR  occurring  within  different  air 
mass  regimes.  The  total  area  compilation  given  in  Table 
6.3  is  useful  for  the  first  three  map  times  because  nearly 
the  entire  grid  area  is  influenced  by  the  same  air 
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mass(see  Fig. 3.1).  However,  at  later  map  times  the  Gulf 
warm  front  moves  northward  and  the  cold  front  moves 
eastward  with  the  result  that  the  domain  becomes 
influenced  by  three  air  masses.  In  order  to  tabulate  grid 
point  totals  representative  of  each  of  these  air  masses 
the  total  grid  was  partitioned  into  three  subregions,  as 
shown  in  Fig. 6. 3.  The  Northeast  subregion  is  defined  as 
those  grid  points  north  of  the  warm  front  that,  developed 
at  OOOOGMT.  The  Southeast  subregion  is  identified  as  that 
part  of  the  grid  south  of  the  warm  front,  and  the  Western 
subregion  encompasses  the  area  behind  the  occluded/cold 
front.  A  simple  count  of  the  number  of  grid  points  where 
precipitation  was  observed  and  latent  heat  release  was 
calculated  in  each  subregion  and  for  each  map  time  is 
presented  in  Table  6.4.  These  statistics  identify  the 
extent  to  which  LHR  was  reasonably  represented  in 
different  air  mass  regions  even  if  the  precise  location  of 
that  LHR  was  not  correct  when  compared  with  the  observed 
precipitation  summaries,  i.e. ,  when  poor  correlations  were 
obtained. 

The  transformation  of  the  observed  precipitation 
pattern  from  scattered  convective  areas  to  a  highly 
organized  synoptic-scale  area  is  readily  shown  by 
comparing  Periods  1  and  2.  In  Period  1 ( 1200GMT- 1 800GMT ) 
the  number  of  points  reporting  observed  precipitation  was 


nearly  the  same  for  the  three  subregions. 
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because  all  three  subregions  are  influenced  by  the  same 
air  mass  type  during  this  period.  As  the  cyclone 
developed  and  moved  eastward,  the  Northeast  subregion 
quickly  gained  predominance,  with  precipitation  occurring 
at  the  majority  of  its  50  grid  points.  The  other  two 
subregions  show  only  modest  increases  in  observed 
precipitation  coverage  after  1800GMT.  The  movement  of  the 
cold  front  into  central  Texas  after  0000GMT  is  shown  by 
the  decrease  in  total  observed  precipitation  points  in  the 
Western  subregion  and  an  increase  in  the  Southeastern 
subregion.  This  table  also  shows  that  the  BLHR 
calculations  did  very  well  in  matching  the  observed 
precipitation  grid-point  totals  in  the  Northeastern  region 
during  Period  2.  However,  during  Period  1  the  BLHR  totals 
are  much  greater  than  the  observed  precipitation  summary 
indicates.  This  overestimation  is  also  noted  in  the  other 
subregions  throughout  the  study  period. 

An  evaluation  of  the  total  number  of  points  where 
latent  heat  release  is  calculated  by  the  modification 
schemes  offers  some  interesting  insights.  In  the 
Northeastern  subregion  the  fo'LHRl  and  MLHR3  schemes  have 
more  points  than  both  the  observed  precipitation  and  BLHR 
at  virtually  every  map  time.  The  MLHR2  modification 
performs  impressively  in  matching  the  number  of  observed 
precipitation  points  in  this  subregion  throughout  the 


study  period. 


88 


89 


7.  SUMMARY  AND  CONCLUSIONS 


This  study  has  investigated  the  utilization  of 
satellite  imagery  to  modify  conventional  latent  heat 
release  estimates.  The  SESAME  I  time  period  of  10-11 
April  1979  provided  a  set  of  enhanced  meteorological  data 
which  could  be  used  to  test  the  well  established  latent 
heat  release  parameterization  technique  described  in 
section  2.3  and  the  use  of  sate  1 1 ite— derived  cloud  cover 
to  modify  these  latent  heating  values. 


An  initial  difficulty  encountered  in  this  study  was  to 
define  a  standard  of  comparison  for  the  various  latent 
heat  release  calculations.  Initial  evaluations  compared 
maps  of  basic  and  modified  LHR  fields  with  hourly  radar 
summaries;  satellite-derived  cloud  cover;  and  gridded, 
tempora 1 1 y— averaged  observed  precipitation  fields.  Then 
the  latter  values  were  used  for  further  evaluation  in 
which  area  averages,  correlation  coefficients,  and  areal 
grid  point  totals  were  examined.  By  combining  these 
approaches  a  more  complete  picture  of  the  LHR  fields  was 
obtained  and  the  impact  of  the  weaknesses  identified  in 
Section  4.1  for  the  individual  comparative  forms  was 
1 essened . 


•>  >  -'■*  A  .*»  *•  . 


During  the  study  period,  a  cyclone  developed  rapidly 
forcing  an  initially  scattered  precipitation  pattern  to 
become  an  organised  synoptic-scale  convective 
precipitation  pattern  after  1800GMT  10  April.  As  shown  in 
the  discussion  in  Chapters  A  and  5,  the  basic  and  modified 
LHP  fields  did  reasonably  well  in  capturing  she  general 
precipitation  pasterns  shown  in  the  radar  and  observed 
precipitation  summaries  in  Period  1 ( 1200GMT— 1800GMT) . 
Yet,  poor  correlation  coefficients  between  the  LHR 
estimates  and  observed  precipitation  values  were  obtained 
for  this  period.  These  low  correlations  were  attributed 
to  the  inability  cf  the  LHR  calculations  to  capture  the 
exact  location  of  the  small-scale  convection  which 
dominated  the  study  area  during  this  early  period.  Still, 
it  is  important  to  note  that  the  modifications  did  improve 
the  LHR  estimates  obtained  from  the  basic  calculations. 
Better  comparisons  between  the  modified  LHR  and  the  radar 
and  observed  precipitation  summaries,  smaller  differences 
between  the  area-averaged  observed  values  and  modified 
LHR.  and  improved  correlation  coefficients  support  this 
c  one  1  us i on . 

As  the  convection  became  more  organized  after  1800GMT, 
there  was  a  large  improvement  in  the  BLHR  calculations. 
This  was  noted  in  the  qualitative  evaluation  but  was 
especially  apparent  in  the  correlation  coefficents. 
Correlations  for  the  BLHR  estimates  increased  from  -0.05 
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at  1800GMT  to  0.46  at  0S00GMT.  Here  again,  the 
modifications  were  able  to  improve  on  the  performance  of 
the  BLHR  calculations.  Furthermore,  correlations  between 
observed  precipitation  and  LHR  calculations  were 
consistently  highest  for  the  MLHR2  scheme. 

Recognizing  that  the  correlation  coefficients  depend  on 
precise  placement  of  calculated  and  observed  values,  a 
second  pattern  comparison  procedure,  which  depends  only  on 
comparability  over  a  general  area,  was  devised.  For  this, 
the  total  grid  area  was  divided  into  three  subregions, 
each  representative  of  a  different  air  mass  region  after 
the  cyclone  and  its  frontal  systems  became  we  1  I  — de ve 1  oped . 
Within  each  subregion  the  number  of  points  with  observed 
precipitation  or  LHR  were  counted  and  compared.  The 
overall  tendency  for  the  MLHRl  and  MLHR3  schemes  to 
increase  the  number  of  points  where  latent  heating  was 
calculated  was  readily  apparent,  as  was  the  decrease  in 
these  LHR  point  totals  resulting  from  the  MLHR2 
modification.  The  increase  by  the  first  and  third 
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to  the  observational  data  forms.  The  use  of  the  235K 
threshold  as  a  discriminator  between  convective  and  stable 
cloud  regions  was  based  on  the  success  that  Arkin(1979) 
had  during  GATE  using  this  temperature  to  correlate 
fractional  cloud  cover  and  6h  rainfall  accumulations.  In 
the  present  qualitative  evaluation  the  MLHR2  modification 
more  consistently  matched  the  precipitation  pattern  and 
the  radar  summaries  than  did  the  basic  and  the  other  two 
modified  LHR  fields.  Furthermore,  the  vertical  cross- 
sections  derived  from  this  modification,  which  had  its 
greatest  impact  on  the  magnitudes  of  the  convective  LHR 
values,  reduced  the  LHR  estimates  to  magnitudes  that  were 
much  more  comparable  to  the  vertical  profiles  of  Fuelberg 
et  a  I .  (1985)  and  the  vertical  cross-sections  of  Robertson 
(1983) . 

The  MLHR2  modification’s  departure  from  the  area- 
averaged  observed  precipitation  totals  were  similar  in 
magnitude  to  the  other  schemes.  However,  the  MLHR2  scheme 
had  the  highest  correlation  coefficients  and  the  best 
overall  percent  improvement  over  the  BLHR  correlation 
coef f icients(45~) .  Finally,  the  MLHR2  modification 
matched  the  number  of  observed  precipitation  points  in  the 
three  subregions  more  closely  than  the  other  calculated 
values  for  nearly  every  map  time  used  in  this  study. 

The  success  of  this  modification  is  tempered  by  several 
key  points.  The  first  is  that  the  MLHR2  scheme 
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The  success  of  this  modification  is  tempered  by  several 
key  points.  The  first  is  that  the  MLHR2  scheme 
occassional  1 y  went  too  far  in  the  elimination  of  LHR 
values  at  grid  points.  Although  it  clearly  matched  the 
radar  summaries  better.  it  underestimated  the  areal 
coverage  of  the  observed  precipitation  fields  at  two  map 
t imes ( 1200GMT  and  1800GMT).  Second,  despite  its  positive 
affect  on  the  BLHR  values,  it  yielded  area—  averaged  LHR 
totals  which  underestimated  the  observed  precipitation 
values,  especially  in  Period  2(2100GMT-0600GMT) ;  and  it 
was  unable  to  produce  good  correlation  coefficients  during 
Period  1 ( 1200GMT— 1800GMT) .  A  third  point  focuses  on  the 
nature  of  the  particular  synoptic-scale  disturbance  under 
study.  This  cyclone  was  a  we  1 1 -deve 1  oped  spring  storm 
with  a  well-established  regime  of  severe  convection  and 
associated  clouds  extending  high  into  the  atmosphere.  In 
cases  containing  significant  precipitation  caused  by 
convective  clouds  with  tops  predominantly  warmer  than 
235K,  this  modification  would  eliminate  much  of  the  latent 
heating . 

The  methods  used  in  this  study  to  derive  and  modify 
latent  heat  release  collectively  represent  a  relatively 
simple  approach  to  a  very  complex  problem.  While  there  is 
little  doubt  that  the  use  of  satellite  imagery  to  improve 
LHR  estimates  shows  merit,  even  more  complex  modifications 
using  the  approach  of  this  study  would  seem  unwarranted. 
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